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ABSTRACT 

Using observations in the COSMOS field, we report an intriguing correlation between the star formation 
activity of massive (~ 10 1L4 M Q ) central galaxies, their stellar masses, and the large-scale (^10 Mpc) envi- 
ronments of their group-mass (~ 10 13 6 M Q ) dark matter halos. Probing the redshift range z = [0.2, 1.0], our 
measurements come from two independent sources: an X-ray detected group catalog and constraints on the 
stellar-to-halo mass relation derived from a combination of clustering and weak lensing statistics. At z = 1, we 
find that the stellar mass in star-forming centrals is a factor of two less than in passive centrals at the same halo 
mass. This implies that the presence or lack of star formation in group-scale centrals cannot be a stochastic pro- 
cess. By z = 0, the offset reverses, probably as a result of the different growth rates of these objects. A similar 
but weaker trend is observed when dividing the sample by morphology rather than star formation. Remarkably, 
we find that star-forming centrals at z ~ 1 live in groups that are significantly more clustered on 10 Mpc scales 
than similar mass groups hosting passive centrals. We discuss this signal in the context of halo assembly and 
recent simulations, suggesting that star-forming centrals prefer halos with higher angular momentum and/or 
formation histories with more recent growth; such halos are known to evolve in denser large-scale environ- 
ments. If confirmed, this would be evidence of an early established link between the assembly history of halos 
on large scales and the future properties of the galaxies that form inside them. 
Subject headings: cosmology: observations — galaxies: groups: general — galaxies: halos 



1. INTRODUCTION 

Understanding the form and evolution of the relationship 
between galaxy stellar mass, galaxy color, and dark matter 
ha lo mass has become a critical topic in galaxy formation. 
In lLeauthaud et all d2012l) (hereafter L12) we combined mea- 
surements of the galaxy stellar mass function (SMF), galaxy 
clustering, and galax y-galaxy lensing in the COSMOS survey 
(IScoville et al.l200"7h to place constraints on the stellar-to-halo 
mass relation (SHMR) at 0.2 < z < 1.0 using a halo occupa- 
tion analysis (HOD). In this paper, we focus on the SHMR for 
massive galaxies, M ga i « 1O 11_1L5 M0, within group-scale ha- 
los, Mi, ~ 10 13 5 M Q , across this same redshift range. Updat- 
ing the L12 results, we now separately constrain the SHMR's 
for galaxies that are actively star-forming and those that are 
passively evolving. We compare these results with a sample 
of central gala xies identified in an X-ray selected COSMOS 
group catalog dGeorge et al.ll201 lb . 
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We define a dark matter halo with as having an over- 
density 200 times the mean cosmic density. All calcula- 
tions assume a flat ACDM cosmology of (f2 m ,(78,f4,«s,/2o) 
0027200807. 0.0438. 0.96.0 0.72) 

2. DATA 

The COSMOS sample that we use for clustering, lensing, 
and SMFs has already been described in detail in L12. HOD 
anlaysis is performed on these measurements. The main dif- 
ference with respect to L12 is that we now divide the sam- 
ple into star-forming ( SF) and passive sub samples using the 
UVJ color-color cuts of lBundv et al.l (120101) . We use the same 
stellar masses as L12. These have b een esti mated us i ng the 
Bayesian code of lBundv etaD (120061) using a lChabrierl (120031) 
initial mass function. In our redshift range, there are 12,573 
passive and 41,682 SF galaxies in the COSMOS sample above 
our completeness limits. 

We also use a COSMOS X-ray selected group catalog 
to select and study central gal axies. Details regarding this 
group catalog can b e found in Finogu enov et al.l (120071) and 
lGeorgeetal T (l2011l) 10 . Halo masses for these groups were 
determined in Leauthaud et al.l d2010l) by calibrating the Lx- 



Mh relation from weak lensing. To ensure a clean sample 
of groups and centrals, we exclude potentially merging sys- 
tems, and groups near maske d regions or with ve ry few mem- 
bers (flag_include=1 in iGeorge et al.l 1201 lT) . This sam- 
ple contains 129 groups out of 211 extended X-ray detec- 
tions. We further remove 18 groups with ambiguous identi- 
fication of a central galaxy, i.e., when the most massive gro up 
galaxy within the NFW-scale radius dNavarro et al.l 1 19971) of 
the halo is not the most massive galaxy within the virial radius 
(mmgg_scale_mstar y mmgg_r200_mstar). At fixed 
redshift, the group catalog constitutes a roughly halo-mass 

10 This group catalog is publicly available and can be found at 
http://irsa.ipac.caltech.edu/data/COSMOS/tables/groups/ 



limited sample of dark matter halos. We divide the data into 
three redshift bins that span z = [0.2, 1.0]. The specific redshift 
bins are the same as in L12 and are shown in Fig. [TJ The mean 
logarithmic halo mass in each redshift bin is 13.47, 13.59, and 
13.75. We note that the mass calibration of iLeauthaud et all 
(2010) assumes a halo mass definition of 200 times the crit- 
ical density. We have converted these values to our fiducial 
halo definition by assuming the NFW density profile with 
a conc entration-mass relation given by iMunoz-Cuartas et all 
d201 lb . We then rescale the masses from the 200-cr itical def- 
inition to the 200-mean (e.g jHu & Kravtsovll2003l) . 

The central galaxies in our sample are well above the com- 
pleteness limit for COSMOS (Fig. 1 in L12), even for passive 
galaxies. We also check for AGN contamination, which we 
will discuss subsequently. 

3. HALO OCCUPATION ANALYSIS 

In ILeauthaud et alj d201 lb . we presented a theoretical 
framework for modeling combined measurements of the SMF, 
galaxy clustering, and galaxy-galaxy lensing. This method is 
a more generalized vers ion of the traditional Ha lo Occupation 
Distribution (see, e.g., ICoorav &~S heth 2002| for a review). 
Our HOD method utilizes these three statistical measures to 
infer the number of galaxies within halos as both a func- 
tion of halo and galaxy mass. In LI 2, we implemented this 
formalism on stellar-mass defined samples within COSMOS. 
Our analysis constrains the halo occupation of both central 
and satellite galaxies, but in this paper we focus exclusively 
on central galaxies within group-scale halos. We constrain a 
SHMR for both passive and SF central galaxies such that the 
total number of central galaxies per halo is unity. This result 
is obtained independent of the SHMR constrained from the 
group catalog. In L12 we assumed that every halo has one 
central galaxy; here we require that the sum of mean occu- 
pation of passive and SF central galaxies is unity. In a com- 
panion paper (J. Tinker et al., in preparation), we present full 
details of our measurements and our model fits. Our results 
focus on the relative clustering of groups at ~ 10 Mpc. Due to 
the small area of COSMOS, the integral constraint can affect 
the clustering of objects at our scale of interest (L12). How- 
ever, it will not alter the relative clustering of two samples in 
the same volume (L12), which is the quantity of interest here. 

4. RESULTS 

4. 1 . Stellar-to-Halo Mass Ratios 

Fig. Q] shows the constraints on the SHMR for passive and 
SF galaxies for each redshift bin. Results from L 1 2 for the full 
stellar-mass limited samples are shown for comparison. At 
z = 0.9, there is a clear difference between the stellar masses 
of SF and passive central galaxies in groups of similar halo 
mass. At M h = 10 13 ' 7 M Q , the difference is 0.4 dex. This is 
qualitatively consis tent with the tren ds seen in AEGIS groups 
at lower halo mass (I Woo et al.ll2.Q12b . At lower redshift, how- 
ever, this difference gradually goes away. In the lowest red- 
shift bin, the SHMR for SF and passive galaxies cross at fixed 
halo mass. 

This evolution is confirmed in the galaxy group sample: 
At high redshift, there is a 3-cr difference in the median cen- 
tral galaxy mass between passive and SF centrals. Errors for 
this quantity are calculated by bootstrap resampling of the 
stellar masses within each subsample. The median galaxy 
masses are also in good agreement with those derived from 
the halo occupation analysis. At lower redshifts, the differ- 
ence in the passive and SF galaxy masses gets monotonically 
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FIG. 1. — Evolution of the stellar-to-halo mass relation for group-scale 
halos. In each panel, the blue and red curves indicate the mean stellar mass 
as a function of halo mass for SF and passive central galaxies, respectively, 
from the HOD analysis. The shaded region around each curve is the 68% 
confidence region on this mean. The green dashed curve shows the result 
for stellar-mass selected samples (no color cut) from L12. In each panel, all 
plot symbols represent results from the group catalog. Small plot symbols 
show the halo and central galaxy masses individual groups; blue stars repre- 
sent star-forming central galaxies, while red circles represent passive central 
galaxies. Objects with X-ray AGN activity are indicated with a black box. 
The larger points with error bars show the median mass of the c entral galax- 
ies in the group catalog, and the uncertainty on that quantity. The green and 
yellow squares show the median values when splitting the sample by disk- 
like (green) and bulge-dominated (yellow) morphologies. We note that these 
results are obtained independently from the HOD results. 



smaller. The median masses do not cross over, as they do in 
the SHMRs; there is a discrepancy between the results for the 
passive subsample in the lowest redshift bin. However, re- 
sults from SDSS demonstrate that this cross over has indee d 
occurred by z = dMandelbaum et al.l2006l:lMore et al.l201 lb . 
The results from the group catalog are qualitatively similar if 
one breaks the catalog up by morphology 1 1 (as shown in Fig- 
ure [T). The galaxies with X-ray AGN activity, either in the 
XMM or Chandra observations are indicated on the plot. The 
low number of such objects, and the (lack of) correlation with 
star formation estimates indicates that AGN contamination is 
not playing any role in the observations. Removing these ob- 
jects from the sample does not shift the medians beyond their 
l-<7 errors. 

At z = [0.22,0.48], the discrepancy in the SHMR values 
and those obtained from the groups for passive galaxies is a 
2.4-er difference based upon creating Monte Carlo samples 
of halos using all elements in the MCMC chain but with the 
same mass distribution as the groups sample. The large-scale 
clustering amplitude of all structure in the low-z bin is be- 
low average (L12), which could drive systematic errors in 
the HOD results. It is also possibl e that the halo mass func- 
tion assumed in the HOD analysis dTinker et a l. 2008a) is not 

11 The 'spheroidal' classification of §3.4.1 in Bundy et al. 2010. 
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the same as the true mass function in that patch of sky, also 
resulting in systematic biases. However, we note that the 
groups catalog is a subset of all X-ray groups within COS- 
MOS, while the SHMRs are derived from statistics on the full 
sample of gal axies. Also, z = measu rements of the SHMR 
usi ng lensing (iMandelbaum et alJ2006h and satellite kinemat- 
ics dMore et al.l201 lb find that the halo masses for massive red 
galaxies are higher than those of massive SF central galaxies, 
following the evolutionary trend seen in the HOD results. 

4.2. Clustering by Central Galaxy Type 

Fig. |2] shows the cross-correlation between groups and all 
galaxies. We split the groups into samples with SF and pas- 
sive centrals, and cross-correlate each set of centrals with all 
galaxies. Both samples are in the z=[0.74, 1.00] redshift bin, 
with a magnitude cut of if&uw = 24. Galaxies brighter than 
this threshold have relia ble photometric red shifts with errors 
~ 0.03 (see figure 2 in iGeorge et al.ll201 ll) . Because most 
of the redshifts of the central group galaxies are known with 
spectroscopic precision (most are sampled within the zCOS- 
MOS survey), we can measure the real-space projected clus- 
tering, which has higher signal-to-noise relative to a simple 
angular cross-correlation. We denote this clustering statistic 
uj(R). Details are given in lPadmanabhan et al.l(l2009l) . Briefly, 
we calculate the projected comoving separation between each 
group-galaxy pair from the angular separation and the redshift 
of the group. We restrict all pairs to lie within a redshift in- 
terval of ±0.09, or 3-cr of the photo-z error. To properly nor- 
malize £j{R) requires detailed information of the photo-z error 
distribution function, but since we are only concerned with the 
relative clustering between two spectroscopic samples cross- 
correlated with the same photometric sample, this step is un- 
necessary. We measure the angular cross-correlation, w(9), 
for the same samples as a cross-check on our results. Errors 
are obtained by bootstrap resampling of the groups and recal- 
culating H)(R) or w(9) for each bootstrap sample. 

A scale of importance is the 1 Mpc scale (comoving), 
roughly the virial radius of the groups. Inside this scale, 
the cross-correlation probes the number of satellite galaxies 
within the groups. Outside this scale, the cross-correlation 
probes the large-scale bias of the groups, which is an indica- 
tor of their environment. In the measurements of Fig. [2] this 
scale marks the bifurcation in the two correlation functions. 
Outside this scale, the passive-central groups have a lower 
large-scale bias, indicating that these halos have formed in 
lower-density environments. Inside this scale, the correlation 
functions for the groups differ at the 1 a level, with the SF- 
centered groups having lower clustering, but the large errors 
prevent meaningful interpretation. 

For each redshift bin, we calculate the bias relative to an 
(arbitrarily normaliz ed) non-linear matte r correlation function 
calculated using the Smith et al. (2003) fitting function. We 
calculate bias using bins at R co > 1 Mpc or 8 > 80 arcsec. 
We estimate the covariance matrix by bootstrap resampling of 
the groups with replacement. We use 200 bootstrap samples. 
Due to the low number of groups, the clustering signal around 
each group can be considered independent. We use the full 
covariance matrix to obtain the bias and its error. The relative 
bias of SF-centered and passive-centered groups is shown as 
a function of redshift in Fig. [3] We show bias measurements 
from both £j(R) and w(8). For the former, a bias measurement 
is not possible at z = 0.36 due to noise in the Q(R) measure- 
ments for both subsamples. At z = 0.88, both measures indi- 
cate that the SF-centered groups have significantly enhanced 
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FIG. 2. — The cross-correlation function of the X-ray groups with all 
galaxies in the defined redshift range. The jc-axis is the comoving projected 
separation between pairs. The y-axis, ui(R), has an arbitrary normalization, 
thus the relative amplitude is the key quantity (see text for details). Black/red 
circles represent groups with passive central galaxies; blue/green squares rep- 
resent groups with star-forming central galaxies. Note that the groups with 
passive centrals are slightly more massive than the groups with star forming 
centrals. 
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FIG. 3. — The relative large-scale bias of SF and passive-centered groups 
(R c0 > 1 Mpc). Circles and squares represent bias obtained from the ra- 
tio of the angular cross-correlation function, w(6), and the real-space cross- 
correlation function, &(R), respectively. Valid bias values cannot be obtained 
for the z = 0.36 u)(R) measurements. The shaded band is the predicted range 
for relative bias of halos separated by angular momentum, with high angular 
momentum halos being more clustered (see the discussion in §5). The dashed 
line indicates no difference in the clustering between the two subsamples. 



clustering. At z = 0.66, the relative bias is above unity, but this 
detection is not significant given the errors. At z = 0.36, the 
angular clustering yields a \ -a detection of elevated clustering 
in the SF-centered groups. 

5. DISCUSSION 

Proposed quenching mech anisms for massive galaxies, 
such as major mergers (i.e.. [Hopkins et"ai1 120081) or AGN 
(i.e., iDe Lucia & Blaizoi l2007h essentially remove galaxies 
from the star-forming sequence and place them in the red se- 
quence. The efficiency of such processes is correlated with 
halo mass. Some research suggests that cooling fl ows may 
lead t o episodic star formation in central galaxies (Liu et all 
12011) . 

The substantial difference in the SHMRs for SF and pas- 
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sive galaxies at the group-scale halo masses at z = 1 has many 
implications. It implies that star formation is not stochas- 
tic in these objects: if massive central galaxies underwent 
periodic episodes of star formation followed by longer-term 
quiescence, the galaxies at fixed halo mass would have the 
same mass regardless of color. The results also imply that 
massive quenched galaxies had far different growth histories 
than those that are forming stars at z = 1. A scenario in 
which galaxies at fixed halo mass grow on a common star- 
forming sequence, with a quenching mechanism that removes 
these galaxies from this sequence, would make passive cen- 
tral galaxies less massive than SF central galaxies. This is the 
opposite of what is observed at z = 1 . To be consistent with 
our observations, passive central galaxies at z = 1 form their 
stars rapidly at high redshift, essentially getting 'ahead of the 
growth curve' relative to central galaxies that are still forming 
stars by z = 1 . At high redshift, central galaxies essentially 
"knew" they would be quenched by z = 1 . 

Fig.|2]indicates that color-selected groups represent special 
subsets of objects at this halo mass scale. The current growth 
rate (indicated by galaxy color) and growth history (probed 
by total stellar mass) of the central galaxy is correlated with 
large-scale environment. A similar effect is seen in dark 
matter halos in N-body simulations, an effect called assem- 
bly bias. For massi ve systems, younger halos exist in mor e 
dense environments dWechsler et al.ll2006t iDalal et al.ll2008l) . 
The environment (and formation history) of massive halos is 
also correlated with angular mom e ntum of dark matter halos 
(IWechsleril200Tt iBett et all 120071: iGao & White! l2007h such 
that high-spin halos are more clustered than low-spin halos. 
This effect goes away below ~ 10 l2 Mp). The shaded r egion 
in Fig. [3]is the numerical result from lBett et al.l (120071) . The 
lower limit is the bias of the top 20% of halos, ranked by 
angular momentum, relative to all halos. The upper limit is 
the bias of the top 20% of halos relative to the lowest 20% 
of halos (see their Fig. 20 12 ). The proper comparison to our 
measurements will lie in between. 

Recent hydrodynamic simulations of galaxy formation in- 
dicate that galaxy morphology is correlated with the angular 
momentum gained from the larger-scale e nvironment around 
the halo at early epochs dSales et al. 201 1). If the angular mo- 
mentum of dark matter and baryons are connected, massive 
halos with central disk galaxies should have enh anced clus- 
tering, in agreement with the results in Fig. |2] ISales et aTl 
d201 lb also find that the disk galaxy masses are lower than 
their spheroidal counterparts. There are, however, substan- 
tive differences between the Sales simulations and our results: 
they show results for less-massive halos at z = 0, rather than 
group-scale systems at z = 1 . They conclude that there is lit- 
tle correlation between morphology and z = halo spin, but 
there is little correlation bet ween halo spin at ear ly and late 
epochs for their halo masses (IVitvitska et al.l2002l) . However, 
the existence of assembly bias implies that more massive ha- 
los retain memory of the angular momentum at the epoch of 
galaxy formation. Further investigation is required at higher 
masses and redshifts. 

However, one need not invoke angular momentum to 
achieve both the rela tive clustering and relative m asses of pas- 
sive and SF centrals. IConrov & Wechslerl (12009b demonstrate 

12 We assume that assembly bias is fixed for halos with the same u(M,z), 
thus we convert their z = results, which are plotted as a function of M, 
to cr(M, z) and interpolate the assembly bias at the values of cr(M, z) for the 
groups samples at each redshift. 



that stellar mass growth peaks at a halo mass of ~ 10 12 M Q , 
weakly dependent on redshift, but the star formation effi- 
ciency at that peak decreases with cosmic time. In this sce- 
nario, central galaxies within late-forming halos would lag 
behind those in early-forming halos, and have enhanced clus- 
tering. This toy model does not explain the morphology de- 
pendence, or the difference in instantaneous SF rates at z = 1, 
but does provide a connection between halo formation history 
and galaxy properties. 

There have been many attemp ts to find a s sembly bias 
in the z = galaxy distribution. iTinker et al.l d2008bl) and 
iTinker et al.l (1201 ll) find no evidence for assembly bias for 
galaxies below the knee in the SMF or luminosity function. 
The assembly biases in the low mass and high mass halo 
populations are driven by disparate physical mechanisms. 
Younger halos form in denser environments at high mass 
through the statistics of Gaussian random fields. At low mass, 
older halos form in denser environments due to tidal forces 
and i nteractions with nearby massive objects dDalal et al.l 
120081) . It is plausible that these two mechanisms may have 
di fferent levels of im pact on galaxy formation. 

IWang et al.l(l2008l) find thatz = group-mass halos in SDSS 
with redder total galaxy content (centrals and satellites com- 
bined) are more clustered than groups with bluer galaxies. It 
is not clear how the clustering signal in z = 1 COSMOS data 
could reverse if the most-clustered halos at one redshift re- 
main t he most clustered at a lower redshift. The IWang et all 
d2008l) detection is mitigated by the lack of independent con- 
straints on the halo mass; in their group-finding algorithm, the 
halo mass is estimated statistically by assuming a 1 : 1 corre- 
spondence bet ween total group stel lar mass and halo mass, 
with no scatter. iBerlind et al.l {2006), using a different group- 
finding algorithm (but the same data set), find the opposite 
signal: groups with bluer central galaxies are the ones that 
are more clustered. Both these methods rely on inferring halo 
mass statistically from the galaxies within them; our X-ray 
detections and lensing masses are more legitimate for detect- 
ing assembly bias. 

From z = 1 to z = 0, the SHMRs evolve quite differently 
depending on star formation activity. By z = 0.36, the mean 
relations have crossed and passive central galaxies live in 
higher mass halos than SF central galaxies at fixed mass. 
Thi s inversion is also consistent with results fr om z = stud- 
ies dMandelbaum et al.ll2006l: iMore et al.ll201 ll) . Star forming 
g alaxies grow by a fac tor of ^2 using the star formation rates 
of lNoeske etaiTd2007|) from z = 0.88 to z = 0.36. Group-mass 
halos also grow by a factor of ~2, thus central galaxies grow 
as fast as their host halos. For quenched galaxies, their growth 
rates are slower than that of their host halos, plausibly causing 
the inversion of the SHMR seen in Fig.Q] 

At z > 1, our results imply that the process that shuts down 
star formation in massive galaxies cannot be explained by a 
stochastic process that is a function of halo mass. Rather, the 
interplay between the dark matter halo and the surrounding 
environment, including the tidal fielde, strongly influences the 
fate of the galaxy forming within it. 
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